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1 Introdução 


Este projeto tinha como grande objetivo introduzir o autor ao ambiente 
de pesquisa, garantindo a ele experiência e conhecimento útil para que 
possa prosseguir pesquisando na área de Teoria de Números ou em áreas 
que se utilizem da Teoria de Números. Os tópicos selecionados para serem 
estudados foram: 

• Propriedades Elementares dos Números Primos 

• Funções Geradoras de Primos 

• Teoremas e Conjecturas ligados a números primos, como o Teorema de 
Dirichlet sobre progressões aritméticas e a Conjectura de Bunyakovsky 

• Criptografia 

A ideia deste projeto de pesquisa nasceu da leitura de Existem funções 
que geram os primos? jl], e ele como um todo gira em torno da busca 
por funções, em especial as polinomiais, capazes de gerar primos com uma 
alta frequência, e por possíveis aplicações destas. Vários textos, que se 
encontram na parte de Referências, foram usados e ainda mais programas 
foram escritos durante o último ano com este propósito, alguns presentes 
neste relatório. 

Durante o período de agosto de 2018 a janeiro de 2019, resumidamente, o 
autor estudou as propriedades elementares dos números primos e as aplicou 
na análise de Polinómios Geradores de Primos. Vários resultados já obtidos 
anteriormente por outros pesquisadores, e que podem ser encontrados nos 
livros de Ribenboim [51 El Ej , foram verificados e alguns foram até mesmo 
atualizados, como o de polinómio quadrático mônico com maior geração 
de primos no intervalo I = [0,1000] D Z do domínio. 

Na segunda metade do período de vigência deste projeto, de fevereiro 
de 2019 a julho de 2019, o foco foi o estudo do RSA e o aperfeiçoamento 
de algoritmos feitos na primeira metade. 

Todo o cronograma proposto foi cumprido. 

2 Materiais e Métodos 

Os materiais utilizados durante a pesquisa foram os livros que se en¬ 
contram nas Referências Bibliográficas e os programas escritos em C pelo 
autor que se encontram abaixo, todos utilizam as bibliotecas stdio.h, st- 
dlib.h, math.h e time.h e foram compilados pelo GCC 8.3. 
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2.1 PrimeLister 


O programa abaixo, PrimeLister, cria um vetor dinâmico já contendo 2 e 
3 como primeiros elementos e vai se atualizando com cada novo número que 
não tem divisores já no vetor. O vetor então é transformado em um arquivo 
list.txt contendo um número MAX de primos a ser usado na maioria dos 
demais programas. 

#define MAX 1000 

int prime(int number, int* v) { 
int raiz=(int)sqrt(number); 
int k=0; 

while(v[k]<=raiz){ 

if (number°/ 0 v [k++] ==0) return 0; 

> 

return 1; 

> 

int main(){ 
int counter=2; 

FILE *list; 

list=fopen("list.txt", "w"); 
int *v; 

if(MAX>=1000000){ 

v=malloc((int)(MAX/10)*sizeof(int)); 
if(v==NULL){ 

printf("malloc has failed\n"); 
system("PAUSE"); 
return 1; 

> 

> 

else{ 

v=malloc(MAX*sizeof(int)); 
if(v==NULL){ 

printf("malloc has failed\n"); 
system("PAUSE"); 
return 1; 

> 

> 

fprintf(list, "2\n3\n"); 
v[0] =2; v[l] =3; 

{int i,j; 

for(i=6, j =2; i<MAX; i+=6){ 
if(prime(i-1, v)){ 

fprintf(list, "°/„d\n",(i-1));counter++; 
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v [j++] =i-l; 

> 

if(prime(i+1, v)){ 

fprintf (list, "7 0 d\n", (i+1)); counter++; 
v [j++] =i+l; 

> 

> 

> 


fclose(list); 
printf ("°/„d",counter); 
return 0; 


2.2 PolyFinderl 

O programa abaixo, PolyFinderl, foi umas das primeiras versões a ser 
usada na busca por polinómios que gerassem bastantes primos no intervalo 
[0,MAX]. Ele ainda tem uma função “prime” que testa a primalidade de 
cada possível primo, feita com base na propriedade de que todo número 
primo maior que 3 é congruente a ±1 mod 6. 

O código é bastante simples: Ele cria um vetor com três elementos, 
representando coeficientes dos polinómios a serem testados, varia os coe¬ 
ficientes de x° e x 1 dentro da coroa circular de raio externo LIM e raio 
interno LIM2 (métrica do supremo), varia o coeficiente de x 2 de 1 a 2, e 
testa a primalidade de cada número gerado pelos polinómios no intervalo 
[0,MAX]. 

O programa não testa alguns polinómios considerados ruins, aqueles 
cujo termo independente é um numero composto, e salva em um arquivo de 
texto PolyList.txt apenas os polinómios que gerarem mais que MINPRIME 
primos. Também é gravado no arquivo o tempo que a busca demorou. 

#define LIM 10001 
#define LIM2 0 
#define MINPRIME 650 
#define MAX 1000 


int prime(int number) { 
if (number<0){ 

number=-number; 

} 

if(number==l){ 
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return 0; 


else if (number<=3){ 
return 1; 

> 

else if (number°/ 0 2==0 I I number°/ 0 3==0){ 
return 0; 

> 

else{ 

unsigned int i,sq; 

sq=sqrt(number); 

for (i=5; i<=sq; i+=6) { 

if (number °/ 0 i == 0 II number 0 /,(i + 2) == 0) 
return 0; 

> 

return 1; 

> 

> 

int main(){ 


clock_t start = clock(); 
int poly[3],i,j,x,y,counter; 

FILE *PolyList; 

PolyList=fopen("PolyList.txt","w"); 


for(poly[2]=1;poly[2]<=2;poly[2]++){ 

for(poly[1]=-LIM;poly [1]<=LIM;poly [1]++){ 

for(poly[0]=-LIM;poly[0]<=LIM;poly[0]+=2){ 

if(poly [0]>-LIM2 && poly [0]<LIM2 && poly[1]>-LIM2 && 
vft poly[l]<LIM2){ 

poly [0]=LIM2; 
continue; 

> 

if(prime(poly[0])==0 && poly[0]!=1)continue; 


f or (counter=0 ,x=0;x<=MAX;x++){ 

y=poly[0] + (x*poly [1]) + (x*x*poly [2]); 
if(prime(y))counter++; 

} 

if(counter>MINPRIME){ 
fprintf(PolyList," ( 0 / 0 d) + (7 0 dx) + (°/ 0 dx''2) : 
m °/od\n" ,poly[0] ,poly [1] ,poly [2] ,counter); 


> 

clock_t stop = clockO; 

float elapsed = (float)(stop - start) / CL0CKS_PER_SEC; 
fprintf(PolyList, "TIME: %f s", elapsed); 
fclose(PolyList); 
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return 0; 

> 

2.3 PolyFinder2 

O programa abaixo, PolyFinder2, é uma versão mais eficiente do pro¬ 
grama acima. O grande diferencial é se utilizar de um arquivo gerado pelo 
PrimeLister para criar um vetor dinâmico que auxiliará nos testes de pri- 
malidade. Deve-se antes saber o tamanho (tamanho lista) desse vetor a 
ser criado. 

Esta versão é em média 35% mais rápida que a anterior de acordo com 
as buscas realizadas. 

#define LIM 10001 
#define LIM2 0 
#define MINPRIME 650 
#define tamanho_lista 350 


int prime(int number, int* v) { 

if(number <0)number=-number ; 
if(number==l)return 0; 

int raiz=sqrt(number); 
for(int i=0;v[i]<=raiz;i++){ 
if (number°/ 0 v [i] ==0) { 
return 0; 

> 

> 

return 1; 

> 

int main(){ 

clock_t start = clock(); 
int poly[3],x,y,counter; 

FILE *PolyList; 

PolyList=fopen("PolyList.txt","w"); 
FILE *list; 

list=fopen("list.txt","r"); 
int* v; 

v=malloc(sizeof(int)*tamanho_lista); 
if(v==NULL){ 

printf("malloc falhou\n"); 
return 1; 

> 
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for(int i=0;i<tamanho_lista;i++){ 
fscanf(list, "°/ 0 d", &v [i]); 

> 

fclose(list); 

for(poly [2]=2;poly[2]<=2;poly[2]++){ 

for(poly [1]=-LIM;poly [1]<=LIM;poly[1]++){ 

for(poly[0]=-LIM;poly [0]<=LIM;poly [0]+=2){ 
if(poly [0]>-LIM2&&poly[0]<LIM2){ 
poly [0]=LIM2; 
continue;} 

if(prime(poly[0], v)==0 && poly[0]!=1)continue; 

for (counter=0 ,x=0;x<=1000;x++){ 

y=poly[0]+(x*poly[l] )+(x*x*poly[2]); 
if(prime(y, v))counter++; 

> 

if (counter>MINPRIME)fprintf (PolyList, " ( 0 /od) + (°/odx) + (°/odx~2) : 
^ °/ 0 d\n", poly[0] , poly[l] , poly [2] , counter); 

> 

> 

> 

free(v); 

clock_t stop = clockO; 

float elapsed = (float)(stop - start) / CL0CKS_PER_SEC; 
fprintf(PolyList, "TIME: %f s", elapsed); 
fclose(PolyList); 
return 0; 


2.4 PrimeCounter 

Há também o PrimeCounter, uma variante mais simplista dos progra¬ 
mas acima, que conta quantos são os primos no conjunto imagem Im = 
(P(0), P(l),..., P(MAX)} de um polinómio quadrático P(n) = ax 2 + bx + 
c. Este não passou por otimizações recentes pois não foi necessário. 

#define MAX 10000 
#define a 1 
#define b -1 
#define c 41 

int main(){ 

int primality=l,counter; 
long long int y,n,i; 

for(n=l, counter=0;n<=MAX;primality=l ,n++){ 
y = (a*n*n + b*n + c); 


6 




if(y%2==0)primality=0; 


for(i=3;primality==l && i<=sqrt(y);i+=2){ 
if (y°/oi==0)primality=0; 

> 

if(primality==l)counter++; 


printf ("There were °/ 0 d prime values for P(n) from n=0 to n=°/ 0 d",counter,MAX); 

return 0; 

> 


3 Resultados 

3.1 Polinómios Geradores de Primos 

Os principais resultados deste projeto, sem dúvidas, foram obtidos du¬ 
rante o recorrente estudo dos Polinómios Geradores de Primos. Com o 
apoio de “Introdução à Teoria dos Números ”[E], que foi extensamente lido 
no primeiro mês de vigência deste projeto, foi possível começar a entender 
mais sobre polinómios como o “Polinómio de Euler”, P(n) = n 2 — n + 41, 
que produz uma sequência de 40 primos distintos para 1 < n < 40 e con¬ 
tinua a produz primos com alta frequência mesmo para valores grandes de 
n. 

Teorema 1. Sep(x ) G Z[x\ e {q,d,r} G Z ; então p(qd+r) = p{r) mod d. 

Analisando os polinómios p q {n) = n 2 — n + q e o Teorema 1, é fácil ver 
que p q (q ) é necessariamente um número composto (ou 1), e que, portanto, 
uma sequência de primos produzida por p q {n ) tem no máximo comprimento 
l = q — 1. 

É natural que se deseje saber quais os polinómios com essa propriedade, 
afinal já se sabia que P 4 i(n), o polinómio de Euler, assim se comportava. 
Com a ajuda de “Prime-Producing Quadratics” [3] foi estudada a caracte¬ 
rização dos polinómios da família p q (n) = n 2 — n+q com máxima produção 
ininterrupta de primos, ou seja, com {p q (l),p q (2), ...,p q {q — 1)} todos pri¬ 
mos. 

Teorema 2. [3] p q {n ) = n 2 — n + q tem máxima produção de primos 
ininterrupta se, e somente se, q G {1, 2,3, 5,11,17,41} 
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O livro “Números Primos, Velhos Mistérios e Novos Recordes” [7] traz 
várias informações sobre a produção de primos dep 4 i(n), como 7r* 2+n+41 (10 7 ) 
2208197. Abaixo definições úteis. 

Definição 1. A função 7r(n) : N —> No retorna o número de primos em 
I n = {1, 2,..., n}. 

Definição 2. A função : Z[x\ xN —X No retorna o número de 

primos no conjunto Im n {\p(x)\) := {|p(l)|, |p(2)|,..., |p(n)|}. 

Com o objetivo de verificar essas informações e adquirir algumas mais 
foi feito o PrimeCounter ( |2.4[ ). Algumas das informações obtidas: 


Tabela 1: Número de primos em I n e em Imj{\'p q {x)\) 


n 

7T (n) 

Kud n ) 

K,U n ) 





10 2 

25 

32 

14 

30 

48 

60 

86 

10 3 

168 

189 

93 

165 

288 

365 

581 

10 4 

1229 

1410 

629 

1208 

2057 

2627 

4149 

10 5 

9592 

10751 

4899 

9086 

15661 

20127 

31985 

10 6 

78498 

88118 

40037 

74058 

128171 

164220 

261081 

10 7 

664579 

745582 

339824 

626458 

1083223 

1388247 

2208197 


Tabela 2: Proporção /n(n) de primos, arredondada nos milésimos 


n 

q = 1 

q = 2 

q = 3 

q = 5 

9= 11 

9 — 17 

q = 41 

10 2 

1.280 

0.040 

0.560 

1.200 

1.920 

2.400 

3.440 

10 3 

1.125 

0.006 

0.554 

0.982 

1.714 

2.173 

3.458 

10 4 

1.147 

0.001 

0.512 

0.983 

1.674 

2.138 

3.376 

10 5 

1.121 

0.000 

0.511 

0.947 

1.633 

2.098 

3.335 

10 6 

1.123 

0.000 

0.510 

0.943 

1.633 

2.092 

3.326 

10 7 

1.122 

0.000 

0.511 

0.943 

1.630 

2.089 

3.323 


Para entender esse comportamento bastam os teoremas [T] e [2j que nos 
dizem que os números gerados por p q (n) não são divisíveis por nenhum 
primo menor que q. 

Com o fim de seguir adiante e descobrir novos formatos de Polinómios 
Geradores de Primos foi feito o PolyFinderl ( |2.2[ ), que encontrou p(n) = 
n 2 — 839n + 6047, um polinómio tal que 7r* 2 _ 839n+6047 (10 3 ) = 693. Segundo 
Ribenboim em seus livros “The Little Book of Bigger Primes” \ 5] e “Die 
Welt der Primzahlen, Geheimnisse und Rekorde ” [6] o recorde de produção 











de primos para um polinómio quadrático mônico no intervalor [ 0 , 1000 ] é 

7r ri 2 -999n+27744l(-*-0 3 ) = 669. 

Na tentativa de encontrar polinómios ainda melhores, o PolyFinderl foi 
aperfeiçoado e tornou-se PolyFinder2 d2.3| ). A partir dele, descobriu-se 

n n 2 -mn+7: 8569(10 3 ) = 704 


. Vale notar que: 


n 2 - 997n + 78569 = (n - —) 2 - 

9 , 839 x o 

n 2 — 839n + 6047 = (n- —) 2 — 


679733 

4 

679733 

4 


Os dois são o mesmo polinómio, deslocado no eixo horizontal. 

Baseando-se nisso foram feitos outros programas similares aos PolyFin- 
ders, que tentavam explorar a simetria vertical das parábolas para en¬ 
contrar melhores geradores de primos. A ideia era encontrar um ótimo 
Polinómio Gerador quadrático para o intervalo [0,500], que tivesse vértice 
próximo à origem, para que então fosse deslocado 500 unidades para a 
direita e se tornasse um ótimo Polinómio Gerador quadrático para o inter¬ 
valo [0,1000]. A busca no intervalo [0,500] é menos custosa que no intervalo 
[0,1000], afinal. Esta busca não encontrou novos bons geradores, mas en¬ 
controu vários já conhecidos, como o n 2 — 997n + 78569. 

Os livros de Ribenboim também falam sobre polinómios recordistas não 
mônicos. A maior produção de primos no intervalo [0,1000] por um po¬ 
linómio se deve a 7r2 X 2_ 158 4 a;+ 9 86 2i = 706. O PolyFinder2 o pôde encontrar 
em suas buscas. 


3.2 Criptografia RSA 

O autor estudou o funcionamento da criptografia RSA com a leitura do 
livro “Números Inteiros e Criptografia RSA” [9j. O conteúdo visto será 
mostrado através de seu uso: 

A chave pública será 

(n, e) = (6028859,5) 

Sabe-se que n = 997 * 6047, números que já apareceram como coefici¬ 
entes de polinómios relevantes neste relatório. Além do mais, 0(n) = 
<£(997)0(6047) = 996 * 6046 = 6021816 e gcd(0(n), e) = gcd(6021816, 5) = 

1 . 
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A tabela de conversão a ser usada será 

Tabela 3: Tabela de Conversão 


A 

B 

C 

D .. 

.. X 

Y 

Z 

10 

~iT 

12 

12 .. 

.. 33 

“34" 

35 


E a mensagem a ser transmitida será IMECC. 

No processo de pré-codificação, IMECC torna-se 1822141212 e então 
18221 41212, pois precisa-se de blocos de informação menores que n = 
6028859. 

Na codificação: 

18221 5 m 1309416 mod 6028859 
41212 5 = 435471 mod 6028859 
Mensagem codificada: 1309416 435471. 

Para decodificar: 

<f>(<f>(n)) = 0(6021816) = 0(2 3 * 3 * 83 * 3023) = 4 * 2 * 82 * 3022 = 1982432 
Disso temos que 

5 1 = 5^ (n))_1 m 5 1982431 = 4817453 mod 6021816 


Para decodificar: 

1309416 4817453 = 18221 5 * 4817453 = 18221 A: *^ n ) +1 ee 18221 mod 6028859 
435471 4817453 3 4 1 2 1 2 5 * 4817453 s 41212**^"^ s‘ 41212 mod 6028859 
graças ao Teorema de Euler. 

O processo de codificação sempre pode ser feito por qualquer um que 
conheça a chave (n, e) e tenha um computador suficiente para fazer as 
contas, mas a decodificação, quando n bem escolhido, só é possível para 
quem conhece a fatoração de n, pois permite calcular 0(n). 

3.3 Demais resultados 

Foram estudados vários teoremas importantes para a Teoria de Números, 
como o Teorema de Dirichlet sobre Progressões Aritméticas, o Teorema de 
Korselt e o Teorema de Green-Tao; e conjecturas como a de Bunyakovsky 
e o décimo problema de Hilbert, todos com o objetivo de conhecê-los e 
formar repertório matemático e histórico. 
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4 Conclusões 


O aluno, durante este projeto, como esperado, adquiriu experiência e 
conhecimento na área de Teoria de Números: Ao estudar os Polinómios 
Geradores de Primos, a Criptografia RSA e demais assuntos relacionados 
ganhou boa base para prosseguir na área e ganhou cultura matemática. 

Vários algoritmos foram desenvolvidos para ajudar nesta pesquisa, cada 
um tendo passado por diversos aprimoramentos, e certamente todos serão 
aprimorados mais diversas vezes quando necessário. Tais algoritmos encon¬ 
traram o polinómio p(n) = n 2 — 997n + 78569, possivelmente o polinómio 
quadrático mônico com maior densidade de primos no intervalo [0,1000] 
até o momento. A busca por polinómios mônicos que gerem mais primos 
em [0,1000] que este p(n) não é fácil, pois teria-se que buscar entre po¬ 
linómios com coeficientes cada vez maiores, que geram números também 
cada vez maiores. Além de aumentar o custo computacional, sabe-se que a 
probabilidade de encontrar primos fica cada vez menor à medida se busca 
entre números maiores. 

5 Perspectivas de continuidade 

O autor tem bastante interesse em continuar a pesquisar sobre Cripto¬ 
grafia, em aprender mais sobre novos métodos criptográficos e poder fazer 
comparações com o RSA. Um novo projeto com objetivo de estudar Crip¬ 
tografia de Curva Elíptica estava previsto para o período de agosto de 2019 
a julho de 2020, foi selecionado em primeira chamada, mas o autor não o 
pôde aceitar pois estará estudando fora do país com uma bolsa provida 
pela UNICAMP no edital n°54/2019 da DERI, logo dará seguimento a 
estes estudos sozinho no momento. 

6 Outras informações 

Foram realizadas pelo autor deste relatório duas palestras essencial¬ 
mente sobre Funções Geradoras de Primos. A primeira, em 14 de setembro 
de 2018 no ciclo de palestras de Curvas Algébricas e temas afins no IMECC, 
sob o nome “Algumas Propriedades da Sequência de Números Primos”. A 
segunda, com título “Polinómios Geradores de Primos”, em 22 de outu¬ 
bro de 2018 pelo DivulgaMat IMECC, um ciclo de palestras voltado para 
alunos da universidade. 
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